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Abstract 



Oh! ^ 

O , A measurement of triple gauge boson couplings is presented, based on W-pair data recorded 

by the OPAL detector at LEP during 1998 at a centre-of-mass energy of 189 GeV with an 

integrated luminosity of 183 pb~^. After combining with our previous measurements at centre- 

X ; of-mass energies of 161-183 GeV we obtain K^=0.97lg:?|], 9l=0.99l1:^om7 and A=-0.110+[5:[li, 

c^ ■ where the errors include both statistical and systematic uncertainties and each coupling is 

determined by setting the other two couplings to their Standard Model values. These results 

are consistent with the Standard Model expectations. 
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1 Introduction 



W-pair production in e^e~ annihilation involves, in addition to the t-channel //-exchange, the 
triple gauge boson vertices WW7 and WWZ which are present in the Standard Model due to 
its non-Abelian nature. Since the start of LEP operation at and above the W-pair threshold 
several measurements have been made of the triple gauge boson couplings (TGCs) involved 
with these vertices in W-pair production [jlj-Q, single W and single photon production 0. 
Limits on TGCs also exist from studies of di-boson production at the Tevatron |0. In this 
paper we update our previous measurements |l|-p to include the data collected during 1998 at 
a centre-of-mass energy of 189 GeV with an integrated luminosity of 183 pb"^. 



According to the most general Lorentz invariant Lagrangian If^^IOl there can be seven in- 
dependent couplings describing each of the WW7 and WWZ vertices. This large parameter 
space can be reduced by requiring the Lagrangian to satisfy electromagnetic gauge invariance 
and charge conjugation as well as parity invariance. The number of parameters reduces to five, 
which can be taken as gl, n^, k^, A^ and A^ [0, §• These parameters are directly related to 
the W electromagnetic and weak properties ||, §]. In the Standard Model gl=hiz=K^=l and 
Az=A'y=0. Precision measurements on the Z° resonance and lower energy data, are consistent 
with the following SU(2)xU(l) relations between the five couplings 0,^, 

Akz = — A/t^tan^ 6*^0 + A(yfJ, 

Az = A-y. 

Here A indicates a deviation of the respective quantity from its Standard Model value and 6*^ 
is the weak mixing angle. These two relations leave only three independent couplings, Ak^, 
Agf and A(=A^=Az) which are not significantly restricted |11, 1^ by existing LEP and SLC Z° 
data. 

Anomalous TGCs give different contributions to different helicity states of the outgoing W- 
bosons. Consequently they affect the angular distributions of the produced W bosons and their 
decay products, as well as the total W-pair cross-section. Ignoring the effects of the finite W 
width and the initial state radiation (ISR), the production and decay of W bosons is completely 
described by five angles. Conventionally these are taken to be [^^,|T^: the W^ production 
polar angleQ 9w', the polar and azimuthal angles, 6** and 0*, of the decay fermion from the W~ 
in the W~ rest frameQ; and the analogous angles for the anti-fermion from the W"*" decay, 6^ 
and 02. However for W-pairs observed in the detector, the experimental accessibility of these 
angles and therefore the sensitivity to the TGCs, is limited and depends strongly on the final 
state produced when the W bosons decay. 

In this study all final states are used, namely the fully leptonic, iuiiui/, the semileptonic, 
qq£z7£, and the fully hadronic, qqqq final states, with branching fractions of 10.6%, 43.9% and 
45.6% respectively. This study is divided into two parts. In the first part the W-pair event rate 



^ The OPAL right-handed coordinate system is defined such that the origin is at the geometric centre of the 
detector, the z-axis is parallel to, and has positive sense, along the e^ beam direction, 9 is the polar angle with 
respect to z and (j) is the azimuthal angle around z. 

^ The axes of the right-handed coordinate system in the W rest-frame are defined such that z is along the 
parent W flight direction and y is in the direction e^ x W where e~ is the electron beam direction and W is 
the parent W flight direction. 



for each of the three final states is analysed in terms of the TGCs. The second part is a study 
of the W-pair event shape for each final state using optimal observables. 

Most of the TGC results of this analysis are obtained for each of the three couplings sep- 
arately, setting the other two couplings to their Standard Model values. These results will be 
presented in tables and log L curvesQ. We also perform two-dimensional and three-dimensional 
fits, where two or all three TGC parameters are allowed to vary in the fits. The results of these 
fits will be presented by contour plots. 

The following section includes a short presentation of the OPAL data and Monte Carlo 
samples and in section ^ the analysis of the W-pair event rate is described. Moving on to the 
event shape analysis, the first step, namely the event reconstruction, is presented in section 
H The TGC extraction using the technique of optimal observables is explained in section |^. 
Systematic errors are described in section ^ and the combined TGC results are presented in 
section ^ Section || summarises the results of this study. 



2 Data and Monte Carlo models 



The data were acquired during 1998 with the OPAL detector which is described in detail 
elsewhere WM . The integrated luminosity, evaluated using small angle Bhabha scattering events 



observed in the silicon tungsten forward calorimeter, is 183.05 ± 0.40 pb ^. The luminosity- 
weighted mean centre-of-mass energy for the data sample is A/i=188.64ib0.04 GeV. 

In the analyses described below, a number of Monte Carlo models are used to provide 
estimates of efficiencies and backgrounds as well as the expected W-pair production and decay 
angular distributions for different TGC values. The majority of the Monte Carlo samples were 
generated at ^/i = 189 GeV with M^^ = 80.33 GeV/c^. All Monte Carlo samples mentioned 



below were processed by the full OPAL simulation program [^ and then subjected to the same 
reconstruction procedure as the data. 

The main Monte Carlo program used in this analysis is the ExcALlBUR |T^ generator. This 
program generates all four-fermion final states using the full set of electroweak diagrams in- 
cluding the W"^W~ production diagrams (clas^ CC03) and other four-fermion graphs, such as 
e"'"e~ — >■ Wel^e, e"'"e~ -^ 7Pe^e~ and e"'"e~ -^ TPTP . Using this program, samples were generated 
with and without anomalous TGCs. Each sample was generated with a different set of TGC 
values. These samples are used to calculate the TGC dependence of the expected event rate 
and angular distributions. We also apply a reweighting technique based on the matrix element 
corresponding to all the contributing four-fermion diagrams. In this way, the angular distri- 
butions for any particular set of anomalous couplings are obtained from Monte Carlo samples 
generated at a limited number of different TGC values. 

For systematic studies we use Standard Model four-fermion samples generated by the 
grc4f |T^ and KORALW [|TB| programs. KORALW uses the same matrix element as grc4f 



^ Throughout this paper, log L denotes negative log-hkehhood. 

^ In this paper, the W pair production diagrams, i.e. i-channel v^ exchange and s-channel 7P /^ exchange, 
are referred to as "CC03" , following the notation of |7t] . 



but has the most complete simulation of ISR out of all the four-fermion Monte Carlo programs 
used in this analysis. The efficiency obtained from the KoRALW sample is used to correct the 
expected event rate. To estimate the hadronisation systematics, Monte Carlo samples were 
produced by grc4f including the CC03 diagrams only, and the fragmentation stage was gener- 
ated separately with either Jetset [|1^] or Herwig ||20|. To study final state interaction effects 



in qqqq events we used Monte Carlo samples with different implementations of Bose-Einstein 
correlations and colour reconnection effects, as detailed in section ^. 

There are other background sources not associated with four-fermion processes. The main 
one, Z°/7 -^ qq, including higher order QCD diagrams, is simulated using Pythia, with 
Herwig used as an alternative to study possible systematic effects. Other background pro- 
cesses involving two fermions in the final state are studied using KoRALZ [^] for e^e~— >-/i"'"/i~, 
e"'"e~— >-r^r~ and e"'"e~— >-i/z7, and Bhwide |^ for e"'"e~— s^e+e". Backgrounds from two-photon 



processes are evaluated using Pythia, Herwig, Phojet [^ and the Vermaseren genera- 



tor |2^. It is assumed that the centre-of-mass energy of 189 GeV is below the threshold for 
Higgs boson production. 



3 Event rate TGC analysis 



The event rate analysis is based on the total W"'"W~ production cross-section measurement ||25[ . 
We use here the same selection with the same efficiency and background estimates. The dif- 
ference is only in the physics interpretation. To analyse the W"'"W~ event rate, all three final 
states are used. Each final state corresponds to a different selection algorithm. These algo- 
rithms are similar to those used at lower centre-of-mass energies (see 0] and references therein) 



and the differences are described in [p5 |. 



The results are summarised in Table |l]. The efficiencies refer to CC03 W+W^ events and 
include cross contamination between the three final states. To calculate the expected num- 
bers of events we use our total integrated luminosity of 183.05 ± 0.40 pb^^ and the total 
W"'"W~ production cross-section value of 16.26 pb corresponding to our centre-of-mass energy 
of 188.64±0.04 GeV and a W mass of Mw=80.42±0.07 GeV/c^ measured at the TevatronQ |2|. 
The cross-section calculation is done using the RacoonWW [^ and YfsWW3 |^ Monte 
Carlo programs which include the most complete 0{a) radiative corrections using the double 
pole approximation method. The results obtained by these two independent programs agree 



within 0.1% p9[. Furthermore the theoretical uncertainty of 0.42% is improved compared with 



the 2% uncertainty of the Gentle [^ semi-analytic program which has been used in our 
previous publications. However, since we have not yet been able to study the full effect of our 
selection cuts with these recent Monte Carlo generators, we prefer not to reduce our assumed 
theoretical uncertainty of 2% which has only a negligible effect on our final results. 

The four-fermion background is split between final states with and without contributions 
from diagrams containing a triple gauge boson vertex. For each of the three final states in our 
signal the TGC-dependent background is calculated from the difference between the accepted 



^ The LEP results for the W mass are not used for the TGC measurement, since they have been obtained 
under the assumption that W pairs are produced according to the Standard Model, whereas W production at 
the Tevatron does not involve the triple gauge vertex. 
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Selected as 


W+W--> eweliyi' 


W+W-^ qq£z7^ 


W+W--> qqqq 


Combined 


Efficiency % 


82.3±1.2 


87.7±0.9 


86.5±0.8 


86.6±0.6 


Signal events 


258±6 


1145±26 


1175±26 


2578±55 


Background events: 
4-fermion, TGC-dep. 
4-fermion, TGC-indep. 

Z77 ^ ff 
Two-photon 


17±2 
7.0±0.3 
4.1±0.3 
0.9±0.9 


27±4 

34±3 

48±6 

3±3 


10±6 

71±10 

245±17 




54±7 

111±9 

297±17 

4±4 


Total background 


29±3 


112±9 


325±21 


466±24 


Total expected 


287±7 


1257±27 


1500±33 


3044±60 


Observed 


276 


1246 


1546 


3068 



Table 1: Expected and observed numbers of events in each W^W final state for an integrated 
luminosity of 183.05 ±0.40 pb"^ at 188.64 dz 0.04 GeV, assuming Mw = 80.42 ±0.07 GeV/c^ and 
Standard Model branching fractions. The efficiency values and the numbers of signal events 
include cross contamination between the three final states. The errors on the combined numbers 
account for correlations between the systematic errors for different final states. 



four-fermion cross-section (including all diagrams) and the accepted CC03 cross-section. The 
background from other final states^ does not depend on the TGCs. 

The errors listed in Table |I] are due to Monte Carlo statistics, luminosity uncertainty of 
0.2%, theoretical uncertainty in the total cross-section (2%), centre-of-mass energy and W 
mass uncertainties (0.04% of the total cross-section each), data/MC differences, tracking losses, 
detector occupancy and fragmentation uncertainties. A detailed description of all these sources 



can be found in 25 . 



The total number of expected events in each final state is consistent with the corresponding 
number of observed events. Therefore, there is no evidence for any significant contribution from 
anomalous couplings. A quantitative study of TGCs from the W-pair event yield is performed 
by comparing the numbers of observed events in each of the three event selection channels with 
the expected numbers which are parametrised as second-order polynomials in the TGCs. This 
parametrisation is based on the linear dependence of the triple gauge vertex Lagrangian on 
the TGCs, corresponding to a second-order polynomial dependence of the cross-section. Conse- 
quently, the expected number of events for each final state also has a second-order polynomial 
dependence on the TGCs. The polynomial coefficients are calculated from the expected number 
of events at different TGC values, as obtained from the corresponding ExcALlBUR Monte Carlo 
samples. For each final state, the corresponding expected cross section is normalised to the 
Standard Model expectation (sum of signal and TGC-dependent background) listed in Table [^. 
The primary reason for this normalisation is the fact that when calculating the numbers in 
Table |I], RacoonWW and YfsWW3 are used for the cross-section and KoRALW for the ef- 
ficiency. This is considered to be more complete than ExcALlBUR. The normalisation factor 
varies between 0.970 {iuii ui' selection) and 0.987 (qqqq selection). 

The probability to observe the measured number of candidates, given the expected value. 



® A small contribution from the I'ei^eff final states which are partly produced by the W fusion diagram is 
included in the TGC-dependent background. 



is calculated using a Poisson distribution. The product of the three probability distributions 
corresponding to the three final states is taken as the event rate likelihood function. The 
systematic uncertainties are incorporated into the TGC fit by allowing the expected numbers 
of signal and background events to vary in the fit and constraining them to have a Gaussian 
distribution around their expected values with their systematic errors taken as the width of the 
distributions. The systematic uncertainties, excluding those on efficiency and background, are 
assumed to be correlated between the three different event selections. 

Data from lower centre-of-mass energies are included, assuming all systematic errors to be 
fully correlated between energies. The corresponding log L curves are used in combination with 
the results of the event shape analysis, which is described in the following sections. The full set 
of results is then presented in Figure |^, where the event rate contributions are shown as dotted 
linesQ. 



4 Event reconstruction 

Starting from the event sample used for the event-rate analysis, a reconstruction is performed 
to extract the maximum possible information on the W production and decay angles, which 
are then used to extract the couplings. Events which cannot be well reconstructed are rejected 
from the sample. 

Three kinematic fits with different sets of requirements are used in the event reconstruction: 

A. requiring conservation of energy and momentum, neglecting ISR; 

B. additionally constraining the reconstructed masses of the two W-bosons to be equal; 

C. additionally constraining each reconstructed W mass to the average value measured at 
the Tevatron, Mw=80.42 GeV/c^ |6|. 



For qqqq events, where all four final state fermions are measurable, fits A, B and C have 4, 5 
and 6 constraints respectively. For qqel/g and cf^V^ events the number of constraints is reduced 
by 3 due to the invisible neutrino. For qqrI7^ events there is at least one additional unobserved 
neutrino from the r decay, but the momentum sum of the track(s) assigned to the r can still 
be used as an approximation to the r fiight direction, relying on its high boost. The r energy 
is left unknown, resulting in a further loss of one constraint. Finally for tP^ Vfj events, where 
none of the leptons is a r, there are two invisible neutrinos. Hence, six constraints are lost and 
requirement C is needed. 

In the following we discuss the reconstruction of each final state separately. 



^ The logL curves are expected to be symmetric around the event rate minimum position which, in most 
cases, is very close to the Standard Model TGC value. However, for the Ak^ parameter, the event rate acquires 
its minimum at positive Ak^ value. The exact location of this minimum differs between tuil v^i qq^^ and qqqq 
events due to contributions from non-CC03 diagrams. Therefore, the overall logL(AK^) function summed over 
the three final states is no longer symmetric, unlike the corresponding functions for the other TGC parameters. 



4.1 Reconstruction of qqivi final states 

Candidate qqez7e and qq/iz7^ events without a well reconstructed lepton track, which were 
included in the sample used for the event rate analysis, are removed from the present sample. 
In this way, each of the events left in the sample has one track that is identified as the most 
likely lepton candidate. For the case of qqrFr events, either one track or a narrow jet consisting 
of three tracks is assigned as the r decay product. 

The electron direction is reconstructed by the tracking detectors and the energy is measured 
in the electromagnetic calorimeters. For muons the momentum is measured using the tracking 
detectors. As explained above, the direction of r candidates can be directly reconstructed 
whilst the energy can only be derived from a kinematic fit. 

The remaining tracks and calorimeter clusters in the event are grouped into two jets using the 



Durham k± algorithm |]3l|. The total energy and momentum of each of the jets are calculated 



with the method described in [^. Kinematic fit A, requiring energy-momentum conservation, 
is then performed on the events. The qqel^e and qq/iz7^ events are accepted if this one-constraint 
fit converges with a probability larger than 0.001. This cut rejects about 2% of the signal events 
and 4% of the background. 

To improve the resolution in the angular observables used for the TGC analysis, kinematic 
fit C is performed on qqel/g and o^fJ^V^ events. The W-mass constraint in this fit allows for the 
finite W-width0. We demand that the kinematic fit converges with a probability larger than 
0.001. For the 4% of events which fail at this point we revert to using the results of fit A. 

For qqrUr events the kinematic fit B is performed and required to converge with a probability 
larger than 0.025. This cut rejects 14% of the signal and 41% of the background. Furthermore, 
this cut suppresses those qqrz/T- events which are correctly identified as belonging to this decay 
channel but where the r decay products are not identified correctly, leading to an incorrect 
estimate of the r flight direction or its charge. The fraction of such events in the qqi"!/,- sample 
is reduced from 18% to 12%. 

Out of the original number of 1246 qq£z7£ candidates used for the event rate analysis there 
are 1075 events left (368 qqez/g, 387 qq/xl^^ and 320 qqTUr) after all these cuts. Assuming 
Standard Model cross-sections for signal and background processes, the remaining background 
fraction is 4.6%. This does not include cross-migration between the three lepton channels. The 
background sources are: four-fermion processes after subtracting the contribution from the 
CC03 diagrams (2.4%), Z77 -^ qq (2.4%), W+W"^ qqqq and iVel'u^' (0.2%) and two-photon 
reactions (0.1%). 

In the reconstruction of the qqiVi events we obtain cos 6'w by summing the kinematically 
fitted four-momenta of the two jets. The decay angles of the leptonically decaying W are 
obtained from the charged lepton four-momentum, after boosting back to the parent W rest 
frame. For the hadronically decaying W we are left with a two-fold ambiguity in assigning the 



* The W mass distribution is treated as a Gaussian in the kinematic fit. However, in order to simulate the 
expected Breit-Wigner form of the W mass spectrum, the variance of the Gaussian is updated at each iteration 
of the kinematic fit in such a way that the probabihties of observing the current fitted W mass are equal whether 
calculated using the Gaussian distribution or using a Breit-Wigner. 



jets to the quark and antiquark. This ambiguity is taken into account in the analyses described 
below. 

In Figure |1| we show the distributions of all the five angles obtained from the combined 
c^iVi event sample, and the expected distributions for Agf = ±0.5 and 0. These expected 
distributions are obtained from fully simulated Monte Carlo event samples generated with 
ExcALlBUR, normalised to the number of events observed in the data. Sensitivity to TGCs 
is observed mainly for cos6'w The contribution of cos6'|, 0^, cos6'*^ and 0*^ to the overall 
sensitivity enters mainly through their correlations with cos^w 



4.2 Reconstruction of qqqq final states 

Using the Durham k± algorithm |^, each selected qqqq candidate is reconstructed as four jets, 
whose energies are corrected for the double counting of charged track momenta and calorimeter 
energies PH]. These four jets can be paired into W-bosons in three possible ways. To improve 



the resolution on the jet four-momentum we perform for each possible jet pairing the five- 
constraint kinematic fit B. We require at least one jet pairing with a successful fit yielding 
a W-mass between 70 and 90 GeV. The W charge is assigned by comparing the sum of the 
charges of the two jets coming from the same W candidate. More details are described in f^. 



The correct di-jet pairing is identified by a likelihood algorithm ^^ which takes into account 
correlations between the input variables. We use as input to our likelihood the W candidate 
mass as obtained from the kinematic fit B, the angles between the two jets coming from each W 
candidate, the difference between the W candidate masses as obtained from the kinematic fit A 
and the charge difference between the two W candidates. According to the W"'"W~ ExcALlBUR 
Monte Carlo, the probability of selecting the correct di-jet combination is about 79%, and the 
probability of correct assignment of the W charge, once the correct pairing has been chosen, is 
about 77%. 

The fraction of events correctly reconstructed can be increased by cutting on the value L 
of the jet pairing likelihood. Its distribution is shown in Figure 0(a). We cut at L > 0.5, which 
minimises the expected total errors on the measured TGCs. 

After all cuts the probability of correct pairing increases to 87% and the probability of 
correct determination of the W charge is about 82%. The selection efficiency is reduced by 
these cuts from 86% to 54%, whereas the background drops from 23% to 12%. In the data, 
889 events survive all cuts out of the 1546 initially selected. The distribution of the W charge 
separation is shown in Figure ^(b). The distributions of the variables cos6'w, cos6'*,^ and 0*.^ 
are shown in Figure 0. 



4.3 Reconstruction of £h'i£ u^i final states 

As explained above, the reconstruction oi tufl vni final states is possible only by using kinematic 
fit C for the case of no r-lepton in the final state. To reject events with r-leptons we use 
the lepton identification algorithm designed for the W branching ratio measurement |]23|. In 
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addition, events with only one reconstructed cone which were accepted for the event rate 
analysis are removed from the sample. As a result of these two cuts, the number of data 
events drops from 276 to 133. The efficiency for iVil v^i events with i, =e,/i drops from 89% to 
74%. The contamination left in the sample from IVgl v^t final states with a r-lepton is 14% and 
the background from other final states is only 0.5%. 

Kinematic fit C has no constraints and turns into solving a quadratic equation. In the ideal 
case of no measurement errors and satisfying the conditions where fit C is valid, namely no 
ISR and narrow W width, one expects to obtain two real solutions corresponding to a two-fold 
ambiguity in the angle set cos^iy, ^\ and 02- There is no ambiguity in the angles Q\ and 
^2 which have a one-to-one correspondence with the lepton momenta. In realistic conditions, 
however, including ISR, finite W width and measurement errors, one might obtain no real 
solution, but a pair of complex conjugate solutions. In fact, only 92 out of 133 events in the 
data have two real solutions, in agreement with the Monte Carlo expectation of 97.8 events. In 
the complex solution case, the nearest real solution is taken by setting the imaginary parts of 
the complex solutions to zero. 

Figure ^ shows the reconstructed cos 6'w distribution. The data distribution agrees with the 
Standard Model expectation. 



5 Event shape TGC analysis 

The angular distributions from all final states are consistent with the Standard Model expecta- 
tion and no evidence is seen for any significant contributions from anomalous couplings. For a 
quantitative study we use the method of optimal observables based on the following dependence 
of the W-pair differential cross-section on the TGCs, 

daiSl.a) = S^^\Q) + ^ai ■ Sf ^(fi) + E"i"i " sif{^), ai = An^, Agl and A, 

where Q are the five phase-space variables, Q= {cos 6^J^, cos 61, (j)l, cos 6 2, 4>2)- For this kind of 
dependence it has been shown [^ that all information contained in the phase-space variables 
fl is retained in the whole set of nine observables Of ^ = SP{Q)/S^^^{Q), Of = SiP{Q)/S^^\n) 
and Oif = of = ^(f (n)/5W(fi) {t,j = 1,2,3). 

The original idea PB|-P^ was to use only the mean values of the first order observables, O . 



However, this causes some ambiguities which can be resolved by inclusion of the mean value of 

— (2) 

the second order observables, O . For one-dimensional fits, where only one TGC parameter, 
ai, is allowed to vary and the others are set to their Standard Model values, we use only the 
corresponding first order, 0\ , and second order, 0\ , observables. In three-dimensional fits, 
where all three TGC parameters are allowed to vary, all nine observables are usedPl. 



7;/ 



The analysis is done separately for the qqi'z/^, qqqq and iuii up final states. The nine 
optimal observables are constructed for each event k with the set of phase-space variables Qk 

^ Trying to use all the nine observables also for the one-dimensional fits does not yield any gain in sensi- 
tivity, since the optimal observables not related to the fitted parameter are not expected to contribute further 
information. 
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using the analytic expression for the CC03 Born differential cross-section to calculate the values 
of S'(°)(r2fc), si (f2fc), S'j-j (fifc) and S^j {Q,k), taking into account the reconstruction ambiguities 
for that particular final state. For the qqqq channel, the uncertainty in the identity of the W"*" 
and the W~ is taken into account by weighting each hypothesis with its probability, determined 
from the Monte Carlo as a function of the charge difference \Qwi — Qw-^ between the two W 
candidates in the event. 

The expected mean values of each observable as a function of the TGCs are calculated 
from four-fermion Monte Carlo events generated according to the Standard Model as well as 
samples generated with various TGC values. The ExcALlBUR matrix element calculation is 
used to reweight the Monte Carlo events to any TGC value required. The contribution of 
background events is also taken into account in the calculation of the expected mean values 
using corresponding Monte Carlo samples. 

An illustrative example is shown in Figure ^ for 0\'^ and 0\^^ using the qq£z7^ final state. 
Figures |^(a) and |^(b) show the distributions of these observables in the data compared with the 
Monte Carlo expectations at the Standard Model and at Ak^=±0.5. The data agree with the 
Standard Model distributions. In Figure |(c) the expected mean values E\Phl^ and i?[Oii^] 
as functions of Ak^ are shown along with the measured mean values O ^ and O . . 

A x^ fit of the measured mean values O to the corresponding expectations £^[(9] (a) is 
performed to extract the couplings from the data. The covariance matrix for the mean values 
is calculated from the Monte Carlo events. In the calculation of the x^ for the q^iT^i channel, 
the contribution of the qqrl'^ final state is calculated separately, in order to account for the 
worse resolution of qqrVr events compared with qqel^e and qq/iz7^ final states. The log L curve 
is obtained by subtracting from the x^ function its minimum value and dividing by two. For 
illustration, the logL function corresponding to the one-dimensional fit to obtain An^ from 
the qqiui data using O^^l and O^l is shown in Figure ||(d). This logL function can be 
decomposed into three contributions coming from O^^l , O^^l and the correlation between 
them. The first two contributions are defined disregarding the correlation. Each of them has 
two minima at equal height. Therefore, using only one of the two observables to calculate the 
log L function results in an ambiguous two-fold solution, whereas using both observables leaves 
only one solution. Figure ^(d) also demonstrates the importance of the correlation between the 
two observables. 

The method is tested with a large number of Monte Carlo subsamples corresponding to the 
Standard Model and to various TGCs. The size of each subsample corresponds to the collected 
luminosity. Summing up the logL curves from the different subsamples corresponding to the 
same coupling gives results which are consistent with the coupling value used in the Monte 
Carlo generation. Also the distribution of the fit results is centred around the expected value, 
but there are some non Gaussian tails. Therefore, it is important to test the reliability of the 
error interval which is estimated in the usual way, to be the region where the logL function 
has a value of no more than 0.5. This is done by calculating the fraction of subsamples where 
the correct value is within the error interval. The error estimate is considered to be reliable 
if the calculated fraction is consistent with 68%, otherwise the logL functions corresponding 
to the subsamples are scaled down by the necessary factor to obtain 68% of the subsamples 
within the error interval. The same scale factor is then applied also on the corresponding log L 
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function of the data. The values obtained for this scale factor are between 0.95 and 1. 

Performing one-dimensional fits to the 189 GeV data, we obtain the results quoted in 
Table 0. The expected statistical errors, which are the average fit errors of the Monte Carlo 
subsamples, are also quoted in Table ||, demonstrating the different sensitivities for the various 
final states and TGC parameters. 

As a cross-check to the mean optimal observable method, two other fit methods are applied 
to the o^iuf events. The first one is the binned maximum likelihood method which has been 
used to analyse our previous data and is fully described in ll|,|l. For the present data we use a 
distribution in all five phase-space variables. This method yields sensitivities comparable with 
those of our main analysis. 

The second cross-check fit for the qq£z7^ results uses an unbinned maximum likelihood. The 
contribution of each data event to the likelihood is the sum over all Monte Carlo events of their 
weights, which are Gaussian functions of the distance in phase-space to that data event. In this 
way, most of the contribution comes from the Monte Carlo events nearest in phase space to the 
data event. The TGC dependence of the likelihood is obtained by using Monte Carlo events 
at different TGC values. This analysis has been applied so far only in a three-dimensional 
phase-space in cos 6'-w, cos 6} and (f)}, and the obtained sensitivities are inferior to those of our 
main method. 

For the qqqq events our fit method is cross-checked with a binned maximum likelihood fit 
to the shape of the cos6'w distribution as done in our previous analysis 0. Since it uses only 
one phase-space variable, this method yields worse sensitivities than our main analysis. 

The results of all these cross-check analyses are compatible with those of our mean value 
optimal observable method. 



6 Study of systematics 

The systematic uncertainties assigned to the event shape TGC results are listed in Table |^. The 
following sources are considered. 

a) Our analysis relies on the predictions of the ExCALlBUR Monte Carlo generator. The 
grc4f and KoRALW generators have a different calculation of the matrix elements and a 
different treatment of ISR. The associated systematic errors are assessed by applying our 
analysis to event samples generated with these Monte Carlo programs. 

b) The Monte Carlo modelling of the efficiency dependence on the event shape is checked 
by comparing data and Monte Carlo distributions of kinematic variables which might be 
affected by the selection cuts, such as the lepton and jet polar angle and energy, as well as 
the angle between the lepton in qq^F^ events and the closest jet. An agreement is found 
between data and Monte Carlo within the statistical errors which are used to estimate the 
systematic uncertainty due to this source. For the iuii ugi final state the event statistics 
are not sufficient to perform this check. The uncertainty for these two charged lepton 
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Ak^ 


Agl 


A 


qqijye 








Without systematics 


n 12+0-28 


-0.029tH?t 


n 1 90+0.075 


Expected stat. uncertainty 


±0.26 


±0.068 


±0.073 


Including systematics 


-0.1ll°;i 


-0.0341°:°^ 


n -|9P;+0.079 


Including 161-183 GeV data 


-0.10l°:i? 


-0.030+H^^ 


n 1 qq+0.071 


Including event rate 


-o.o7+°:?i 


-0.028t°:[Ji 


n 10^+0.068 


qqqq 








Without systematics 


0.36^^6 


0.081°:!^ 


0.051°:?^ 


Expected stat. uncertainty 


±0.57 


±0.15 


±0.19 


Including systematics 


o.sitn 


'J-U0_o.l5 


o.o6l°:i° 


Including 172-183 GeV data 


0.43t°i^ 


0.121^:11 


o.i7lg:iS 


Including event rate 


0.85tni 


o.iGtrii 


0.30l°i^ 


iuel uei 








Without systematics 


-o.iotSf 


-o.o4l°:i 


-0.19 ±0.15 


Expected stat. uncertainty 


±0.50 


±0.27 


±0.17 


Including systematics 


-0.10+Sf 


-o.o4l°:i 


-0.19 ±0.15 


Including 183 GeV data 


n 1 7+0-53 


-o.i7l°:i 


-0.20 ±0.14 


Including event rate 


-o.o8+°:i; 


-0.091°!^ 


n 1 7+0-13 


All final states 








Including systematics 


-0.03;°:?^ 


-0.018+H^I 


n 1 1 0+0.O66 

— U.ii»_o.064 


Including previous data 


-o.oo+°:i 


-o.oo8tHi 


Q 1 1 q + 0.059 


Overall combined results 








(including event rate) 


-o.o3t°:?° 


-o.oo9t°:r7 


-U.iiU_o.o55 


95% C.L. limits 


[-0.32, 0.45] 


[-0.12, 0.11] 


[-0.22, 0.01] 


3-D fit results 


0.02+°:?° 


n 1 on+o-O'"^ 

U.iZU_o.o83 


-0.l90t'oZ 



Table 2: Measured values of the TGC parameters for each final state and the results after 
combining all final states. For the individual final states we also list the results before including 
the systematic errors and the expected statistical errors. We also list the results after combining 
with our previous data and after combining with the event rate information from all centre-of- 
mass energies. For the overall combined data we list also the corresponding 95% C.L. limits. 
In the last row we present the results of a 3-dimensional fit to the overall combined data where 
all three TGC parameters are allowed to vary simultaneously. 
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events is taken as twice the uncertainty of the qq£z7^ events with one charged lepton. 
This conservative estimate is based on the assumption that the qqii^i uncertainty in the 
efficiency is entirely related to the lepton identification. 

Uncertainties in estimating the background are determined by varying both its shape 
and normalisation. The shape predicted by Pythia of the Z'^/7 — > qq background which 
is the main one for the qqiui and qqqq final states is replaced by that predicted by 
Herwig. The normalisation is varied within the background uncertainty as determined 
in the cross-section analysis p5[ . The contributions from two-photon processes are much 



smaller and a conservative estimate of their uncertainties is obtained by removing them 
from the analysis. 

d) The uncertainty due to jet and lepton reconstruction is investigated by considering the 
effects of smearing and shifting jet parameters in Monte Carlo events. The resolutions 
in jet energy, cos 6* and are varied by 10%, 5% and 3% respectively. Jet and lepton 
energies are shifted by 0.5% and jet cos^ values are shifted by 0.0003. The sizes of 
these variations have been obtained from extensive studies of back-to-back jets in Z° 
events collected during calibration runs. The same reconstruction uncertainties are also 
assumed for the r jets in qqrVr events. The momentum resolution of electrons and muons 
and their charge misassignment probability are modified by varying the resolution in q/pr 
by 10%. Here q and pr are the lepton charge and transverse momentum with respect to 
the beam. 

e) The uncertainty in the beam energy of 20 MeV and the effect of our Monte Carlo samples 
being produced at a centre-of-mass energy 360 MeV higher than the data are investigated 
by applying our analysis to Monte Carlo samples generated at different centre-of-mass 
energies then those used for reference samples. The effect of the W-mass uncertainty is 
assessed in a similar way using Monte Carlo samples generated with different W-masses. 

f) Possible dependences on fragmentation models are studied by comparing the results of 
the analysis when applied to a Monte Carlo sample generated with grc4f using either 
Jetset or Herwig for the fragmentation phase. The main effect is on the qqqq final state, 
where Herwig predicts a higher probability of correct jet pairing and correct W charge 
assignment than Jetset. This effect is partially due to the average charged multiplicity 
for u,d,s,c quarks being lower in our tuned version of Herwig than in Jetset. Since the 
data are in agreement with Jetset, we weight the Herwig Monte Carlo events so as to 
reproduce the same mean charged multiplicity as Jetset. The effect of this weighting is 
a reduction of the systematic error by about 25-30% for all couplings. 

g) Bose-Einstein correlations (BEC) in qqqq events might affect the measured W charge 
distribution. We investigate this effect with a Monte Carlo program simulating BEC via 



re-shuffling of final state momenta |39] 



h) The jet reconstruction and the measured W charge distribution might also be affected 
by colour reconnection. This effect is investigated with several Monte Carlo samples 



generated according to the Ariadne |^ and Sjostrand-Khoze |^ models. The largest 
effect is observed in the Ariadne 2 and Ariadne 3 models. However, these models have 
been disfavoured by studies of three jet events in LEP Z^ data |H2[. For this reason the 



current Ariadne implementation of colour reconnection is not used to assign a systematic 
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uncertainty. Of the remaining models considered, the Sjostrand-Khoze model I produces 
the largest bias. 



Source 


qqii^e 


qqqq 


iuei u^i 


Any 


^9l 


A 


Ak^ 


^9l 


A 


Ak^ 


^9l 


A 


a) MC generator 


0.020 


0.001 


0.007 


0.05 


0.01 


0.01 


0.13 


0.03 


0.02 


b) Efficiency 


0.038 


0.010 


0.012 


0.03 


0.01 


0.02 


0.06 


0.02 


0.02 


c) Background 


0.030 


0.002 


0.007 


0.04 


0.02 


0.03 


0.02 


0.02 


0.02 


d) Reconstruction 


0.040 


0.005 


0.014 


0.01 


0.01 


- 


0.01 


— 


— 


e) -Ebeam and Mw 


0.021 


0.006 


0.014 


0.05 


0.01 


0.03 


0.01 


0.01 


0.01 


f) Fragmentation 


0.008 


0.012 


0.005 


0.17 


0.06 


0.08 


- 


- 


- 


g) BEG 


- 


- 


- 


0.04 


0.02 


0.03 


- 


- 


- 


h) Colour reconnection 


- 


- 


- 


0.01 


0.01 


0.03 


- 


- 


- 


Combined 


0.070 


0.018 


0.026 


0.20 


0.07 


0.10 


0.14 


0.05 


0.04 



Table 3: Contributions to the systematic uncertainties in the determination of the TGC 
parameters from the 189 GeV event shape data for different final states. 

The dominant contributions to the systematic errors in different final states come from 
different sources, and therefore we neglect any possible correlations between the systematic 
errors from different final states. 

The systematic uncertainties from all sources are combined in quadrature. The values in 
Table |] are not directly used in the analysis. Instead we use the systematic uncertainties in 
the expected mean optimal observables. Correlations between the systematic errors of different 
observables are also taken into account and the full systematic covariance matrix is added to 
the statistical covariance matrix. The combined covariance matrix is used in the x^ fits to 
obtain the TGCs. In this way, the TGC errors obtained from the fit are already the combined 
(statistical and systematic) errors. The results of these fits are listed in Table |^. 



7 Combined TGC results 



The logL functions obtained from the TGC event shape analyses including the systematic 
errors are combined with the curves obtained from our previous analysis at lower centre-of- 
mass energies |I|-§[. The systematic errors from the results of the different centre-of-mass 
energies are assumed to be uncorrelated. This assumption does not affect the final results, 
since the small contribution of the systematic uncertainties to the total errors of the previous 
data (see Tables 8, 10 and 12 in [^) becomes completely negligible after combining with the 
present results. 

Table ^ lists the event shape results for each final state after combining with the previous 
data. The corresponding logL functions are added to those obtained from the event rate. The 
correlation between the systematic errors of the event rate and event shape, mainly due to the 
uncertainty in the background level, is neglected, as it affects the results by less than 1% of 
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their statistical errors. The resulting log L curves are shown in Figure ^ and the corresponding 
TGC results are listed in Table ^ 

Finally the results of the three final states are combined. This is done by first combining 
separately the three logL functions from the event shapes and those from the event rate. In the 
next step the event shape and event rate logL functions are combined into the overall result. 
In combining the event shape information from the three channels, the correlations between 
the systematic errors are neglected as already explained in section ^. On the other hand, the 
correlations in the event rate data are not neglected, as already mentioned in section ^ The 
log L curves obtained after combining the three final states are shown in Figure |^, separately 
for the event rate and event shape information and their sum. In Table ^ we list the TGC 
results obtained from the combination over the three final states of the event shape data and 
those after including the event rate information. The overall combined results are shown also 
in terms of 95% CL. intervals. 

To study correlations between the three TGC parameters we also extract a log L as a function 
of all three variables, Ak^, Ag\ and A. Figure H shows the 95% CL. contour plots obtained from 
two-dimensional fits, where the third parameter is fixed at its Standard Model value. We also 
perform a three-dimensional fit, where all three couplings are allowed to vary simultaneously. 
These fit results are listed in the last row of Table |^ and the corresponding two-dimensional 
projections are plotted as dashed contour lines in Figure |[ As can be seen, the allowed range for 
each parameter is extended when the constraints on the other two parameters are removed. The 
different minima locations for the two- and three-dimensional fits in some of these plots, as well 
as their non-elliptic shape, emerge from the presence of local minima in the three-dimensional 
log L function [|T^ . 



8 Summary 

Using a sample of 3068 W^W^ candidates collected at LEP at a centre-of-mass energy of 
189 GeV, we measure the TGC parameters. For this measurement we use the event rates and 
event shapes of W pairs decaying into all final states. After combining our results with those 
obtained from the 161-183 GeV data we obtain: 

^ _ n 07+0.20 

91 = 0.9911°:°^?, 

A = -0.1101°:°!, 

where each parameter is determined when the other two parameters are set to their Standard 
Model values. These results supersede those from our previous publications ||l]-Q] . They are all 
consistent with the Standard Model. 
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Figure 1: Distributions of the kinematic variables cos^w> cos^|, cos^*,^, 0| and 0*,^, as obtained 
from the qqii^e events. The sohd points represent the data. The histograms show the expecta- 
tion of the Standard Model and when Agl = ±0.5. The shaded histogram shows the non-qqiu^ 
background. Notes: 1. In the case of W'^^iu decays the value of(j)} is shifted by vr in order to 
overlay W^ and W^ distributions in the same plot. 2. The jet with < 0*^ < vr is chosen as 
the quark (antiquark) jet from the decay of the W~ (W~^) according to our convention. 
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Figure 2: For the qqqq channel, (a) distribution of the output of the jet-pairing hkelihood 
corresponding to the most likely combination for the data (points) and for the Excalibur four- 
fermion Monte Carlo (histogram) . The hatched area represents the contribution of wrong 
pairing (assuming all four quark final states compatible with W pair production are coming 
from such events), the double hatched area represents the contribution of the TGC-independent 
four-fermion background and the dark area represents the contribution of the Z°/7 -^ qq back- 
ground. The arrow indicates the cut value, (b) Distribution of the charge separation between 
the two W candidates for the data (points) and for the four-fermion Monte Carlo (histogram). 
The hatched area shows the contribution of correct pairing and incorrect W charge, and the 
double hatched area shows the contribution of wrong pairing. The dark area represents the 
contribution of the Z^/'j -^ qq and TGC-independent four-fermion background. 
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Figure 3: Distribution of the kinematic variables cos^w> cos^*^, and 0*,^. for qqqq events. The 
sohd points represent the data. The Monte Carlo predictions for the Standard Model and 
for Agl=+0.5, —0.5 are shown respectively as solid, dashed and dotted lines. The dark area 
represents the contribution of the Z'^/7 -^ qq background and the hatched area represents the 
contribution of the TGC-independent four-ferniion background. The jet with <0*,t< tt is 
chosen as the quark (antiquark) jet from the decay of the W~ (W'^) candidate according to our 
convention. 
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Figure 4: Observed cos^w distribution in the tP^v^i analysis. All events enter with a total 
weight oi one, but in the case of events with two ambiguous solutions for cos^W; each solution 
enters with a weight of 0. 5. The histograms show the expectations of the Standard Model and 
when A = ±1. 
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Figure 5: qqiVf event distributions of the (a) first and (b) second order optimal observables 
corresponding to Ak^. The histograms show the expectation of the Standard Model and when 
Ak.^ = ±0.5. The expected mean values of these observables are plotted in (c) (lower and upper 
plot respectively) as solid curves. The measured mean values and their statistical errors are 
indicated by the shaded horizontal bands. In figure (d) the contributions of the first and second 
order optimal observable to the log L function, disregarding the correlation between them, are 
plotted by the dashed and dotted lines respectively. The dash-dotted line is the contribution 
of the correlation between the two observables and the solid line is the overall log L function 
which is the sum of the three contributions. 
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Figure 6: Negative log-likelihood curves obtained from the different final states: q(\ivi (dashed 
lines), qqqq (dotted lines) and tuii up (dash-dotted lines). Each curve is obtained by combining 
the results from the event shapes and the event rates and setting the other two TGC parameters 
to their Standard Model values. Systematic errors are included. The solid line is obtained by 
combining the three final states. 
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Figure 7: Negative log-likelihood curves obtained using different sources of information on the 
TGCs. The curves for each TGC parameter are obtained setting the other two parameters 
to their Standard Model values. The dashed lines are obtained from the event shapes and 
the dotted lines from the event rates. The systematic errors are included. The solid lines are 
obtained by combining both sources of information. 
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Figure 8: The two-dimensional correlation contours corresponding to AlogL = 2.995 for differ- 
ent pairs of TGC parameters, as obtained in two and three parameter fits. The value 2.995 is 
chosen because in two dimensions it corresponds to a 95% C.L. interval. For the two parameter 
fit (solid curve), the third parameter is fixed at zero. In the case of the three parameter fit, the 
third parameter is not restricted and the dashed curve is a projection onto the two-dimensional 
plane of the envelope of the three dimensional AlogL = 2.995 surface. The solid (open) points 
indicate the best two (three) parameter fit values and the stars indicate the Standard Model 
expectations. These results are obtained from all event rate data as well as event shapes of all 
qq£z7^ data, 183-189 GeV qqqq data and 189 GeV iViJ u^' data. 
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